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ABSTRACT 

We present an initial study of the mass and evolutionary state of a massive and distant cluster, 

RCS2 J232727.6-020437. This cluster, at z=0.6986, is the richest cluster discovered in the RCS2 
project. The mass measurements presented in this paper are derived from all possible mass proxies: 

X-ray measurements, weak-lensing shear, strong lensing, Sunyaev ZePdovich effect decrement, the 
velocity distribution of cluster member galaxies, and galaxy richness. While each of these observables 
probe the mass of the cluster at a different radius, they all indicate that RCS2 J232727.6-020437 is 
among the most massive clusters at this redshift, with an estimated mass of M200 ~ 3 x lO^/i^pM©. In 
this paper, we demonstrate that the various observables are all reasonably consistent with each other 
to within their uncertainties. RCS2 J232727.6-020437 appears to be well relaxed - with circular and 
concentric X-ray isophotes, with a cool core, and no indication of significant substructure in extensive 
galaxy velocity data. 

Subject headings: galaxies: clusters: individual (RCS2 J232727.6-020437) 


1. INTRODUCTION 

High redshift clusters have been successfully identified 
in dedicated surveys working with a range of cluster- 
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selecting techniques and wavelengths. These include 
cluster discoveries in deep X-ray observations (e.g. 

& Luppino|il994[ IR.osati et al? 1 '"' 


Romer et al. 


Gioia 

2000 


Ebeling et al.|2001[ 

Rosati et al. 2004; 

Mullis et al. 2005; 

Stanford et al. 2006 

Rosati et al. 2009 

; optical+near in- 

frared imaging (e.g.,[Gladders & Yeej2005 Stanford et al. 

2005 

Brodwin et al. 2006; |Eiston et ai. 2006; Wittman 

et al. 

2006; |Eisenhardt et al. 2008; Muzzin et al. 2009; 

Wilson et ai. 2009; 

Papovich et al.||2010[ Brodwin et al. 

2011 

Santos et ai. 

2011) Gettings et ai. 2012; Stanford 


and detection ot the the Sunyaev Zel’dovich (SZ) "ef¬ 


fect flStaniszewski et al. |2009 Vandcrlinde et ah] 2010 


K a 

R 


Williamson et al.|2011| |Marriag eTt al. 2011; PianckW 


laboration et al. 120131 [H assclficl Tet al. 1120131 IReichardt 

et al.||2013| Brodwin et ai. 2014 Bleem et al.|2015|). 

Nevertheless, despite this extensive effort, these sur¬ 

veys resulted in a modest number of high redshift (z > 
0.5) and massive (M > 1O 15 M 0 ) galaxy clusters. This 
relative paucity of distant massive clusters is a reflection 
both of the challenges inherent in d etecting s uch c luster s 
and of their intrinsic rarity (e.g., Crocce et al. 2010). 
Such clusters are the earliest and largest collapsed halos; 
the observed density of distant massive clusters is thus 
exqui sitely sensitive t o several cosmological parameters 
(e.g., Eke et al. 1996) and indeed the presence of a sin¬ 
gle cluster m prior cluster s urveys has been used to limit 
cosmological models (e.g., Bahcall & Fan 1998). Such 
clusters also offer, at least in principle, the opportunity 
to test for non-gaussianity on clust er scales if the cosmol¬ 
ogy is otherwise constrained (e.g., Sartoris et al. 2010). 

At 0.6 < z < 1.0, the most massive ga laxy clus- 
ters known to date are CL J1226+3332 (Maughan 
et al.| |2004|) at z = 0.89 with mass of M 2 oo = 
1.38 ± 0.20 x lO 15 /i7nM 0 (|Jee & Tyson||2009i), ACT- 
CL J0102-4915 at z = 0.87 and with M 2 oo = 2.16 ± 
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0.32 x 10 15 /i 70 1 Mq (Menanteau et al. 20121, and 
MACS0744.8+3927 at z = 0.6 98 with M\< 1.5Mpc) = 
2.05 ± 0.57 x lO 15 h 7o 1 M 0 (Applegate et al. 2014). 
Recently discovered z > 1.0 galaxy clusters appear 
to have more moderate masses, e.g., SPT-CL J2106- 
5844 (z = 1.14, M 9 nn = 1.27 ± 0.21 x IO^/i^Mq; 
Foley et al.|[20ll| , SPT-CL J2 040-4451 (z = 1.48, 
M 2 qo = 5.8 ± 1.4 x 10 14 /i 70 1 M q ; Bayliss et al. 2014bI 
and IDCS J 1426.5+3508 (z = 1.7 5, M 200 = 4.3 ± 1.1 x 
Brodwin et al. 2012). 
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Massive clusters at any redshift are amenable to de¬ 
tailed study with a density of data that less massive sys¬ 
tems do not pr esen t. The X-ray luminosity of clusters 
scales as M 180 ( Pratt et al.||2009 ), the SZ decrement as 
M 1 - 66 (Bonamente et al. 2008), the weak-lcnsing shear 

approximately as A/^qq, and the galaxy richness in a 
fixed metric aperture (and hence the available number 
of cluster galaxy targets for spectroscopic and dynami- 
cal studies within a giv en field of view) scales as M 0 6 
(Yee & Ellingson 2003) at these masses. Similarly it is 
expected that the most massive clusters d ominate the 
cross-sec tion for cluster-scale strong lensing (Hennawi et 
aL|2007). Thus the most massive clusters offer a wealth 


of potentially well-measured observables which can be 
used, for example, to study the correspondance between 
different mass proxies; such study is critical to the suc¬ 
cess of surveys which aim to use the redshift evolution of 
the cluster mass function as a cosmological probe. 

We present here detailed observations of a single mas¬ 
sive cluster selected from the Second Red -Sequence Clus¬ 
ter Survey (RCS2; [cTlbank et al. |2011f. This cluster, 
RCS2 J232727.6-020437 (hereafter RCS2327), was se¬ 
lected from RCS2 in an early and partial cluster cat¬ 
alog. Its optical properties indicated that it is a very 
massive cluster, and justified an extensive followup cam¬ 
paign with ground-based and space-based observatories 
at all wavelengths, from X-ray to radio. Since its dis¬ 
covery, some of the properties of R CS2327 have been 
reported on in the literature. Gralla et al. (2011) first 
measured its mass from Sunyaev Zetdovich array obser¬ 
vations and its Einstein radius from strong lens model¬ 
ing. RCS2327 was rediscovered as the highest signifi¬ 
cance cluster in the Atacama Cosmology Telescope sur¬ 
vey (Hasselfield et al. 2013), and Menanteau et al. (2012) 
also report on mass estimates from archival optical and 
X-ray observations. 

Although the discovery publication of RCS2327 has 
been delayed, it was advertised in the past decade in 
various oral presentations and conferences - in order to 
motivate more extensive followup effort by the commu¬ 
nity. Indeed deeper and more detailed observations have 
been conducted since, and will be the basis of future pub¬ 
lications. This paper presents mass estimates from the 
initial survey and early multi-wavelength followup obser¬ 
vations of RCS2327, which collectively indicate that it is 
an unusually massive high-redshift cluster of galaxies. 

This paper is organized as follows. The appearance of 
the cluster in the RCS2 data and catalogs is discussed in 
§ [2] We describe the various datasets and corresponding 
analyses (richness and galaxy photometry, dynamics, X- 
ray, SZ decrement, weak- and strong-lensing) in detail in 
We discuss the implications of these observations in 


§ [4] and conclude in § [5} 

Throughout the paper we use the conventional nota¬ 
tion M 2 oo (A/ 500 , 4 / 2500 ) to denote the enclosed mass 
within a radius R200 (-R500, -R2500), where the overden¬ 
sity is 200 (500, 2500) the critical matter density at the 
cluster redshift. Unless otherwise stated, we used the 


WMA P 5-year cosmology parameters (Komatsu et al. 
20091, with Ua = 0.73, tt m — 0.27, and H 0 = 70 km s -1 


Mpc -1 . In this cosmology, 1" corresponds to 7.24 kpc at 
the cluster redshift, 2 = 0.6986. Magnitudes are reported 
in the AB system. 


2. THE SECOND RED-SEQUENCE CLUSTER SURVEY AND 
THE DISCOVERY OF RCS2327 

The Second Red-Sequence Cluster Survey (RCS2) is 
an imaging program executed using the Megacam facil¬ 


ity at CFHT. RCS2 is described in full in Gilbank et 
al. (2011). In short, images have been acquired in the 


g, r, and z filters, with integration times of 4, 8 , and 6 
minutes respectively, and all with sub-arcsecond seeing 
conditions via observations in queue mode. The RCS2 
data are approximately 1-2 magn itudes deeper th an the 
Sloan Digital Sky Survey imaging (York et al. 2000), with 
a 5-cr point-source limiting magnitudes of 24.4, 24.3, and 
22.8 mag in g , r, and z respectively. The RCS2 sur¬ 
vey data comprise 785 unique pointings of the nominally 
1 square degree CFHT Megacam camera; the surveyed 
area is somewhat less than 700 square degrees once data 
masking and pointing overlaps are accounted for. The 
cluster and group catalog from RCS2 extend t oz~ 1.1, 
constructed u sing the techniques described in |Gladders| 
& Yee (20051. The 5 -band imaging improves the over- 
ail performance at lower redshifts (compared to RCS1; 
Gladders & Yee 2005), and makes the survey more adept 
at detecting strong lensing clusters, since lensed sources 
tend to have blue colors. 

RCS2327 was discovered in 2005 in an early and par¬ 
tial version of the RCS2 cluster catalog. An examination 
of the RCS2 survey images made it clear that it was an 
unusually massive object. A color image of RCS2327 
is shown in Figure [l] The original RCS2 imaging data 
clearly showed at least one strongly lensed arc, and the 
indicated cluster photometric redshift was £ ~ 0.7. A 
plot of the detection significance versus photometric red¬ 
shift for clusters from the RCS2 cluster catalog is shown 
in Figure [2] The RCS2 imaging data are fairly uniform, 
and so at a given redshift the detection significance is a 
meaningful quantity that is not strongly affected by data 
quality from region to region of the survey. At high red- 
shifts RCS2327 is the most significant cluster detected. 
Furthermore, a cluster of a given richness and compact¬ 
ness (both of which influence detection significance) will 
be detected as a more significant object at lower red- 
shifts; the fact that RCS2327 is detected with a signifi¬ 
cance as great as any lower redshift clusters implies that 
it is likely the most massive cluster in this sample. Even 
from these basic data and considering the volume probed 
it is apparent that RCS2327 is a remarkably massive clus¬ 
ter, worthy of significant followup. 

The cluster is located at R.A.=23:27:27.61 (J2000) and 
Decl.= —02:04:37.2 (J2000); this is the position of the 
brightest cluster galaxy (BCG) and is coin ciden t with 
the center of the cluster X-ray emission (see §|3.3|below). 
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Fig. 1.— A composite gri color image of RCS2327 from imaging from the LDSS-3 instrument on Magellan (see § |,3. 1 [ . The field of 
view is 2x2 proper Mpc at the redshift of the cluster, centered on the BCG. RCS2327 is obviously demarcated by the abundance of red 
early-type galaxies in the center of the image. I mmed iately south of the foreground bright galaxy to the north-northeast is an obvious 
strongly lensed merging double image (see Section |3.8[) which was also readily apparent in the original RCS2 survey data. 


3. FOLLOWUP OBSERVATIONS AND MASS ESTIMATES 

In this section we describe the multi-wavelength fol¬ 
lowup observations of RCS2327. Based on these observa¬ 
tions, we are able to estimate the mass of RCS2327 from 
richness, galaxy dynamics, X-ray, Sunyaev Zel’dovich ef¬ 
fect, weak lensing and strong lensing. The different mass 
proxies naturally measure either spherical mass or a pro¬ 
jected mass along the line of sight (usually referred to 
as cylindrical mass or aperture mass). Moreover, each 
mass proxy is sensitive to mass at a different radial scale: 
strong lensing measures the projected mass density at 
the innermost parts of the cluster, typically ~100-500 
kpc, and is insensitive to the mass distribution in the 


outskirts; SZ decrement and X-ray measure the mass at 
larger radii (typically R2500) and lack the resolution at 
the center of the cluster; weak lensing reconstructs the 
projected mass density out to R200, with poor resolution 
at the center as well. Dynamical mass (from the veloc¬ 
ity distribution of cluster galaxies) is used to estimate 
the virial mass. We note that these mass proxies are not 
always independent, and rely on scaling relations and as¬ 
sumptions. In the following subsections, we describe the 
data and our analysis to derive the cluster mass from 
each mass proxy. In § [4] we compare the masses derived 
from the different mass proxies. 
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Photometric Redshift 


Fig. 2.— The detection significance versus redshift of RCS2327 
compared to the distribution of these properties for clusters from 
the RCS2 cluster catalog. RCS2327 is indicated as a heavy symbol 
and is an obvious outlier; it is as significant as any cluster detected 
at lower redshifts. 

3.1. Deep Multi-color Imaging, Galaxy Distribution and 

Richness 



18 19 20 21 22 23 24 

i Magnitude 


In addition to the RCS2 survey imaging data, avail¬ 
able imaging data on RCS2327 includes images from the 
LDSS-3 imaging spectrograph on the 6.5m Clay tele¬ 
scope, taken during a run in September 2005. Total inte¬ 
gration times were 16, 12, and 10 minutes in the g,r and 
i filters respectively. The point-spread-function width 
at half maximum in the final stacked images is 0760 (i), 
0765 (r), and 0780 (g) with some image elongation due to 
wind shake present principally in the bluest band. These 
data cover a circular field of view 8’ in diameter, centered 
on the BCG. Figure |T| is constructed from these data. 

The RCS2 data are best suited to measurements of 
cluster richness, as they are well calibrated and natu¬ 
rally include excellent background data, and are readily 
connected to the cosmological context and calibration of 
the m ass-richness relation provided by the RCS 1 pro¬ 


gram (|Gladders & Yee||2005| |Gladders et al.||2007|). The 
multi-band LDSS-3 images are deeper, with better see¬ 
ing than the RCS2 images, and we use these data for 
computing detailed photometric properties - principally 
color-magnitude diagrams. For this we focus our analysis 
below on the r and i observations, since this filter pair 
has the best image quality, and almost perfectly straddles 
the 4000A break at the cluster redshift. 

Figure [3] shows an r — i color-magnitude diagram of 
all galaxies at projected radii less than 1 Mpc from the 
cluster center. The red sequence of early-type cluster 
members is obvious, emphasizing the extraordinary rich¬ 
ness of this cluster in comparison to most other c lus¬ 
ters in the literature at a similar reds hift (e.g.JDe Lu¬ 
cia et al.||2007| |Gladders & Yee||2005|). Figure |d| shows 
only galaxies for which a spectroscopic redshift is avail¬ 
able, plotted by SED type. As expected, galaxies which 
are both cluster members and have early-type spectra 
are almost all red-sequence members. From the spec¬ 
troscopically confirmed early-type cluster ga laxies, with 
simpl e iterated 3-a clipping (e.g., as in Gladders et al. 
19981, we fit a linear red sequence relation, given by 
(r — i) = —0.053 x (* — i*) + 1.436. We take the char- 


Fig. 3.— The color-magnitude diagram for a 1 Mpc radius field 
of view centered on RCS2327, in the r and i filters. Galaxies have 
been divided into three radials bins of equal radius; galaxies in the 
outermost bin are plotted as small pluses, and those in the central 
bin as heavy squares. Only galaxies with photometric uncertainties 
in both filters of less than 0.2 magnitudes are shown. 



Fig. 4.— The color-magnitude diagram for galaxies with spectro¬ 
scopic redshifts, regardless of position, in the r and i filters. Spec¬ 
troscopic early-type cluster members are shown as red diamonds; 
other cluster members are shown as green squares. Non-members 
are shown as small squares. Only galaxies with uncertainties in 
both filters of less than 0.2 magnitudes are sho wn. T he best fitting 
red-sequence model, derived as described in § |3.l[ is overplotted 
with the characteristic magnitude i* indicated by the heavy black 
square. 
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acteristic magnitude for cluster galaxies in RCS2327 as 
i* = 21.3, consi stent to within 0.05 mag nitud es with the 


models in both Gladders & Yee (20051 and Koester et 
al. (20071. These models include a correction tor passive 


evolution. The measured scatter of early-type galaxies 
about the best fit red sequence is less than 0.05, con- 


of redshifts ( 

Hao et al. 

2009 


Mei et al.| 

2009 

). The 

best fitting model is indicated 

in Figure 141 Tfiese data 


demonstrate that RCS2327 appears as expected for a 
well-formed high-redsliift cluster, albeit an extraordinar¬ 
ily rich example. 

We derive a total richness for RCS2327 of 
B 9 cT= 3271T488 Mpc^ 8 , and a corresponding red- 
sequence richness of -B oc _bs=2590±413 Mpc^ 8 (see 
Gladde rs fe Yee| |2005[ for a detailed explanation of 
the B gc parameter). Calibrations relevant to the 
measurement of this richness have been taken from the 
RCS1 survey, which also uses the z'-band as the reddest 
survey filter. A direct co mparison of the total ri chness 
to the scaling relations in Yee & Ellingson (2003) nom¬ 
inally corresponds to a mass of M 2 oo = 4.6I^J x 10 15 
/lyjMg with a significant uncertainty , given the known 


scatter in richness as a mass proxy ( Gladders et al.|2007 


Rozo et al. 2009), and the lack of direct calibration o: 
the richness-mass relation at the redshift of RCS2327. 
Furthermore, the relevan t richness to use in comp arison 
to the scaling relation in Yee fe Ellingson (|2003 l is not 
obviou s; though the richness values in | Yee fe Ellingson| 
( 2003]) are for all galaxies, the small blue fraction in that 
sample and the significant observed evoluti on in the 
general cluster blue fraction (Loh et al.||2008) from the 
redshift of that sample (mean z=0.32) To the redshift of 
RCS2327 suggests that the (less evolving) red-sequence 
richness may be a more appropriate measure. With that 
in mind we note that the mass corresponding to B gc ns 

is Mom = 3.2~}"n 


- 0.8 


x 1O 15 /iy O 1 M 0 . 


Given the limitations of 


this analysis however, we do not use a richness-derived 
mass extensively in the analysis in § [4j but simply note 
here that RCS2327 is remarkably rich. 

3.2. Optical Spectroscopy 

Spectroscopic observation of galaxies in the field of 
RCS2327 has been conducted using the Magellan tele¬ 
scopes. Data were acquired in both normal and nod-and- 
shuffle modes using LDSS-3, during runs in August and 
November, 2006, and a total of 3 masks with the GISMO 
instrument in June 2008. RCS2327 was also observed us¬ 
ing the GMOS instrument on the Gemini South telescope 
in queue mode in semester 2007B, yielding redshifts of 
potential lensed so urce s; the Gemini data are discussed 
in more detail in § |3.8| below. 

All Magellan spectra have been reduced using standard 
techniques, as implemented in the COSMOS pipclint]^] 
The bulk of the LDSS-3 observations (apart from a sin¬ 
gle early mask, which established the cluster redshift at 
z ~ 0.70) were acquired using a 6000A-7000A band- 
limiting filter; this allows for a high density of slits, at the 
expense of a significant redshift failure rate (specifically, 
[O II]AA3727 is undetectable outside of 0.61 < z < 0.87, 
and the Ca H and K lines are undetectable outside 

21 http: / / obs.carnegiescience.edu/Code/cosmos 


of 0.53 < z < 0.76). The GISMO observations were 
conducted using a band limiting filter covering 5700A- 
9800A. 

A total of 353 robust redshifts were measured from 
these data. Most are unique, with overlap between ob¬ 
servations with different instruments or runs amounting 
to a few galaxies per mask. From six galaxies in common 
beween the LDSS-3 and GISMO data the mean difference 
in redshifts is measured to be 135 km sec -1 , and the un¬ 
certainty within observations using a single instrument is 
measured to be less than 100 km sec -1 . Neither of these 
uncertainties is significant in the analysis below. Red¬ 
shifts were measured using a combination of cross corre¬ 
lation and line measurement techniques, and cross cor¬ 
relation measurements of absorption systems were only 
retained if (at minimum) the H and K lines were individ¬ 
ually visible. Apart from possible mis-interpreted single 
emission line redshifts in the LDSS-3 spectra, the mea¬ 
sured redshifts are robust. Each spectrum was also clas¬ 
sified as either an emission or absorption type, with post- 
starburst (showing strong Balmer lines) or AGN features 
also noted when present. 

A histogram of the galaxy velocities around the mean 
cluster redshift of 0.6986±0.0005 is shown in Figure [6] 
The measured velocity dispersion for 195 cluster mem¬ 
bers is 1563±95 km sec -1 with uncertainties measured 
from a bootstrap analysis. The cluster is well separated 
from other structures. The velocity dispersion using 
only the 110 galaxies with early type spectra is 1398±99 
km sec -1 , and similarly using all other cluster members 
we derive 1757±139 km sec -1 - a factor of 1.27±0.14 
larger. These differences are as expected and in line with 
that observed for relaxed X-ray selected clusters at lower 
redshifts, where the typical ratio in velocit y dispersio n 
of blue to red cluster members is 1.31±0.13 (Carlberg et 


al. 


1997). 


’he velocity distribution of cluster members versus 
projected radius, in Figure [ 7 ] shows several trends also 
consistent with that expected for a relaxed cluster. The 
velocity dispersion is a declining function of cluster¬ 
centric radius, an effect most apparent in the early- 
type galaxies. In non-overlapping radial bins of 0.5 
Mpc in radius, and at mean radii of 0.27, 0.71 and 1.27 
Mpc, we find velocity dispersions of 1626±127 km sec -1 , 
1268±147 km sec -1 , and 1034±201 km sec -1 , respec¬ 
tively. The radial distribution of emission line members 
relative to absorption line cluster members is also as ex¬ 
pected, with proportionately more actively star-forming 
systems found at large radii. A clear interpretation of 
this result is difficult given the complexity of the sam¬ 
pling from multiple masks from multiple instruments 
with differing fields of view, and the weighting of slit 
assignments toward photometric red-sequence members; 
the data shown in Figure [ 7 ] are at least consistent with 
expectations. Finally, a KS test of the velocity distri¬ 
bution shows at best marginal evidence for velocity sub¬ 
structure, with the velocity distribution inconsistent with 
a normal distribution at a modest 1.3 sigma using all 
galaxies. Using only early type members, there is not 
even marginal evidence for velocity substructure. 

We also note the presence of a secondary structure sep¬ 
arated from the main cluster by 5700 km sec -1 . This 
structure is dominated by emission line galaxies, has a 
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Fig. 5.— The distribution of red-sequence galaxy light (red contours) smoothed with a 150 kpc FWHM gaussian, and the full band X-ray 
light (blue contours) similarly smoothed. Blue diamonds mark X-ray point sources; the correspondance between some of these sources 
to optical counterparts validates the astrometric matching between these data. Individual galaxies are indicated by squares, with symbol 
size proportional to brightness; solid red squares have colors consistent with the cluster red-sequence and are the basis for the overplo tted 
contours, and solid yellow squares show spectroscopically confirmed members of the secondary structure noted in Figures [6| and [7] and § |3.2| 
The inset shows the central l'xl'region with the same data, with the strong lensing main halo mass distribution also indicated (green 
contours). The contour levels have been chosen to highlight the core position and outer shape of the light distributions. 


velocity dispersion of 400 km sec -1 , and is located to the 
edge of the spectroscopic field of view, as can been seen 
in Figures [5j [6j and [ 7 ] It is not significant for any of the 
analyses below. 


3.2.1. Dynamical Mass Estimates from Velocity Dispersion 

The observed velocity dispersion is converted to mass 
through t he vir ial scalin g relation derived from simula¬ 
tions ( Evrard et al.|20 08|. This relationship is described 


as 


-^200 — 


10 


( &obs \ 

h(z) V CT15 ) 


(1) 


where a 0 b s is the observed 1-d velocity dispersion of the 
cluster, (t 15 is the velocity dispersion normalized for a 
10 15 M^ cluster, and a is the slope of the scaling rela¬ 


tion. 


Evrard et al. (2008) find the best fit parameters 


for a multitude of cosmologies and velocity dispersions 
measured from dark matter particles to be <Ji 5 = 1082.9 
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Fig. 6. — The velocity distribution of galaxies about the mean 
cluster redshift for all galaxies (top), absorption line only galaxies 
(middle), and all emission line galaxies (bottom). Velocity disper¬ 
sions in the restframe in km sec -1 are as indicated. 



Cluster-centric Radius (Mpc) 


Fig. 7.— Redshifts versus cluster-centric radius for absorption- 
line-only galaxies (red diamonds) and all other galaxies (green 
squares) for all galaxies within A z = 0.05 of the RCS2327. The 
field of view of the GISMO instrument is approximately 0.7 Mpc in 
these coordinates; the denser sampling this provides in the cluster 
core, combined with the details of slit placements in the LDSS-3 
instrument yields a complex slit weighting with radius which is 
responsible for the apparent deficit of galaxies at ~1 Mpc. 


km s 1 and a = 2.975. Using this scaling relation with 
the measured velocity dispersion for all galaxies, the re¬ 
sulting mass is M 2 oo = 2.97^0^5 x IO^/i^Mq. The 
uncertainties in mass assume a 13% total uncertainty 
in velocity dispersion. This included both the statis¬ 
tical uncertainty, which is small for the large number 
of galaxies observed in this sample, and the system¬ 
atic uncertainty, which includes line-of-sight effects, clus¬ 
ter shape/triaxiality, and foreground/background con¬ 
tamination (Gifford et al. 2013, Saro et al. 2013). 


Saro et al. (2013) re-fit the scaling relation to a semi- 
analytic galaxy catalog for the Millennium Simulation 
and find the parameters to be ays = 938.0 km s -1 and 
a = 2.91. The resulting mass using these parameters is 
M 2nn = 3.08tJ;$ x 10 15 WM„. 
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3.2.2. Dynamical Mass Estimates from the Caustic Method 

The distribution of radial velocities of cluster galaxies 
as a function of cluster-centric radius can be used to es¬ 
timate its mass using the caustic technique (Diaferio & 
Geller 1997, Gifford et al. 2013, Gifford & Miller 2013). 
This method relies on the expectation that cluster galax¬ 
ies that have not escaped the potential well of the clus¬ 
ter halo occupy a well-defined region in a radius-velocity 
phase space confined by the escape velocity from that 
poten t ial, v^Ar). We follow the techniques outlined in 
Gifford et al. (2013), and refer to that publication for a 
full description oFthe methods applied here. 

Figure [8] shows the radius-velocity space of ~ 250 
galaxies. We fit an iso-density contour to the data to 
find v esc as indicated by the velocity edge in phase-space 
density. The enclosed mass can be derived as 


r r 200 


GM'200 — 


fFp(r)A 2 (r)dr, 


( 2 ) 


where A 2 (r) is the square of the line-of-sight escape ve¬ 
locity, and J>(r) is a function of the potential, density, 
and velocity anisotropy, corrected for projection effects. 
We apply the common convention of assuming that Tq 
is constant. Physically, this parameter depends on the 
unknown concentration and velocity anisotropy profile 
of the cluster. Disagreement on the constant value of 
Ta that results in unbiased mass estimates on average 
persi sts in t he literature w ith values ranging from 0.5 to 
0.7. Gifford et al. (2013) find that = 0.65 results 
in mass estimates with less than 4% mass bias for sev¬ 
eral semi-analytic catalogs available for the Millennium 
Simulation, and we adopt this value for this study. 

We derive a dynamical mass of M 2 qo = 2.94+J'°g x 
10 15 /i/ 0 1 Mq. The uncertainty in the derived mass using 
the caustics technique depends on the number of galaxies 
used; from caustic mass analysis of t he M illenium Simu¬ 
lation semi-analytic galaxy catalogs, Gifford et al. (2013) 
find that for N ga i > 150 the scatter is < 30% with a bias 
of < 4%. 

We note that the velocity dispersion of galaxies iden¬ 
tified as possible members by the caustic technique is 
1586 ± 58 km s -1 , in agreement with the estimates in 


3.3. Chandra Observations, X-ray properties and Mass 

Estimates 

RCS2327 was observed on two separate occasions with 
the Chandra X-ray Observatory. A first 25 ks observation 
was carried out on 2007 August 12 (Cycle 8 Proposal 
08801039; PI: Gladders) using the ACIS-S array. The 
early analysis suggested that the cluster was massive, X- 
ray regular, and possibly hosting a cool core. These hints 
justified the need of deeper data obtained in 2011 with 
ACIS-I (Cycle 13 Proposal 13800830; PI: Hicks). 

The two deep Cycle 13 ACIS-I pointings (150 ks total) 
sample a more extended field around the cluster, and 
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Fig. 8.— The projected radius vs line-of-sight velocity phase 
space for galaxies identified as members by the caustic algorithm 
(red circles) and non-members (open circles). The estimated escape 
velocity surfaces (or caustics) are represented by solid red lines 
which are symmetric with respect to the cluster velocity. 


result in a higher signal-to-noise ratio than the Cycle 
8 observation. Since the background is better under¬ 
stood than that of the ACIS-S configuration, combined 
with the very small increase in signal-to-noise ratio that 
would be gained by combining both datasets, we chose 
to analyse the 2011 data separately and not co-add the 
two epochs. 

The X-ray data reduction follows Martino et al. (2014) 
and Bartalucci et al. (2014), with minor modifications. 
We use both the count statistics and the control on sys- 
tematics offered by a multi-component mod eling of the 
background noise (see|Bartalucci et al.|2014[ for details). 
Filtering the hard and soft event light curve reduces 
the total exposure time by 10% to ~ 130 ks. We bin 
the photon events in sky coordinates with a fixed an¬ 
gular resolution of T/4 and a variable energy resolution 
that matches the detector response. The effective ex¬ 
posure time and the estimated background noise level 
are similarly binned. Following the Chandra Calibration 
database CALDB 4.6.1, when computing the effective ex¬ 
posure we take into account the spatially variable mirror 
effective area, quantum efficiency of the detector, CCD 
gaps, bad pixels, and a correction for the motion of the 
telescope. Our background noise model includes Galac¬ 
tic foreground, cosmic X-ray background, and false de¬ 
tections due to cosmic ray induced particles. For the par¬ 
ticle backgro und spectrum, w e use t he analytical model 
proposed by |Bartalucci et al.| (|2014[) . The amplitude of 
all the other components was determined from the data 
outside the region of the field of view covered by the 
target. We derive the temperatu re map follow ing tech¬ 
niques described in Bourdin & Mazzotta ([2008]). Figure|9| 
shows the X-ray flux isophotes overplotted on the false- 
color temperature map. The cluster has a regular X-ray 
morphology and d oes n ot show sig nific ant subst ructure 
or X-ray cavities (Hlavacek-Larrondo et al. 2014). 


3.3.1. X-ray surface brightness, temperature, and metallicity 

profiles 

We extract the surface brightness profile (Figure 10 1 ) 
from an effective exposure and background-corrected soft 


band ([0.5-2.5] keV) image, after excluding point sources. 
The profile averages the surface brightness in concentric 
annulli centered on the maximum of a wavelet-filtered 
image of the cl ust er. The temperature and metallicity 
profiles (Figure [l0|3,c) were calculated in five radial bins 
out to ~ 840 kpc, each containing at least 2000 counts in 
the [0.7-5] keV band. The measurements of temperature 
and metallicity assume redshifted and absorbed emission 
spectra modeled with the Astrophysical Plasma Emission 
Code (APEC, Smith et al. 2001) adop ting the element 
abundances of|Grevesse & Sauval| (|1998D and neutral hy- 
drogen absorption c ross sections of Balucinska-Church fe| 


McCammon (1992). The spectra, modified by the effec- 
tive exposure and background, are also convolved with a 
function of the redistribution of the photon energies by 
the detector. The assumed column density value is fixed 
at 4.73 x 10 20 cm 2 from measurements obtained near our 
target by th e Leiden/Argentine/B onn (LAB) Survey of 
galactic HI (Kalberla et al. 2005). The redshift is fixed 
to z — 0.6986. 

We find that the metallicity increases from ~ 0.2 
solar at i ?2500 to ~ 0.6 solar at the cluster core. 
RCS2327 shows a temperature gradient towards the cen¬ 
ter of the cluster, indicating a significant cool core. The 
temperature in the estimated [0.15 — 1] x R 500 region 
(roughly out to IMpc, see § 3.3.2) is Tx = 13.9l 2 'g 


keV. We also estimate the cooling times as a function 
of cluster-centric ra dius (follo wing the p rescription de¬ 
scribed in Hlavacek-Larrondo et al. 12013), and find that 


the cooling time profile decreases mildly towards the cen¬ 
ter from ~ 40 Gyr at 4 00 kpc , to ~ 4 Gyr at the core. 
From the gas density (§ 3.3.2) and temperature profile, 
we find that the central entropy is 46 ± 28 keV cm -2 . 


3.3.2. X-ray Mass Estimate. 

To measure the X -ray mass we follow the f orwa rd pro- 


et al. 


cedur e descr ibed in Meneghetti et al. 
(21)12). In short, analytic mocle 


(2010 


s are 


and Rasia 
itted to the 


projected surface density and temperature profiles and 
subsequently analytically de-projected. The 3D infor- 
mati on are then folded into the hydrostatic mass equa- 

~l5b. 

d i 


tion (Vikhlinin et al. 2006 


parametrized via a modifie' 


The surface brightness is 
/3 —model with a power-law 


trend in the center and a steepening behavior in the out¬ 
skirts, plus a second /3—model to describe the core: 


n p n e (r ) = rig 


(r/r c y 


fc 02 


(l + r 2 /r2)3/3-<*/2 (l + r 2/ r 2 2 )3/3 2 ’ 


( 3 ) 


where n e and n p are the electron and proton densities, 
respectively. We allow all parameters to vary. 

We model the temperature profile with a simple power- 
law: 

TMr) = Mr/r t y a . (4) 


This profile is then projected along the line of sight using 
the formula of the spectroscopic-like temperature: 


_ / WT 3B dV 
f WdV 


( 5 ) 


where W = (n p n e ) / (T^ 5 ). All the best fit parameters 
are determined using a x 2 minimization technique ap¬ 
plied to the models and the data. Finally, the 3D density 
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and temperature profiles are used to estimate the total 
gravitational mass through the equation of hydrostatic 
equilibrium (HSE; jSarazin 19881: 


M(r) = -3.68x10 u T{r)r 


„2 (d log p g , d\ogT\^_ lM 


V dr 


dr 


i 70 JU Oi 


( 6 ) 

where the numerical factor includes the gravity con¬ 
stant, proton mass, and the mean molecular weight, 
H = 0.5954. 

To estimate the uncertainties, we produce 500 realiza¬ 
tions of the surface brightness and temperature profiles 
assuming a Poisson distribution for the total counts in 
each annulus and a Gaussian distribution for the pro¬ 
jected temperature. The fitting procedure described 
above is repeated each time. The derived mass pro¬ 
files are related to the original HSE mass profile via 
the resulting value of the least-mean-square formula 
(= £[M rea i iza tion - -^original] 2 , where the sum extends 
to all radial bins). We consider the 68% of the profiles 
(340 in number) with the smallest associated value and, 
finally, for each radial bin, we consider the maximum 
and minimum values of the selected profiles. The re¬ 
sulting ma ss profile and its uncertainties are plotted in 
Figure [l0] h The HSE radius and mass at overdensity 
A =2500 are -R 2500 = ^^- 33^0 kpc, and M 2500 = 
3.2jtg;3 x 10 14 /i 7 O 1 Mq, respectively, and the gas mass 
within this radius is M gaSi 2500 = 4.41°'® x IO^/i^Mq. 

On average, HSE masses are expecte d to be b iased low 
by 10-15% as evid ent from simulations (Rasia et al. 2006 


Nagai et al. 20071 Battaglia et al.|2013|) and observations 

(|Mahdavi el al. 2008, 2013). This offset is smaller than 
our statistical uncertainty. 

The value of R 2500 is within the region probed by the 
observation, and thus its measurement is conservative 
and robust. However, the lower overdensities A = 500 
and A = 200 are not within the observed region, and 
we therefore need to extrapolate. For that purpose, we 
follow two different approaches. 

1 : NFW-mass extrapolation — We fit the 500 realizations 
with the Navarr o-Frenk-White (NFW, |Navarro et al. 
19951 (19961 (1997]) formula: 


Mnfw = 47rM 0 r s 


log(l +r/r s ) - 


r/r s 


1 + {r/r s )_ 


, (7) 


where r s is the scale radius and Mq the normalization of 
the mass profile. The fitting is carried out only within 
the observed radial region. The results lead to R 500 = 
\.\ 5 to J 2 phd Mpc and R 2 00 = 1.78t^/ifg 1 Mpc. The 
errors represent the minimum and maximum values of 
the 68% of the analytic expressions that are the clos¬ 
est to the NFW fit of the original mass profile. The 
resulting extrapolated masses are M 500 = I.IIqJ x 
IO^/j^Mq and M 200 = 1.8±J;| x IO^/i^Mq. 

2: M — Yx relation. — To derive R 500 we also apply the 
iterative method bas ed on the M — Yy(= M yas x Tx) 
relation proposed by |Kravtsov et al.| (|2006|) ■ We start 
with an initial guess for the radius (we consider twice 
the value of the measured i? 25 oo)- We evaluate the gas 
mass at that radius from the surface brightness profile 
and compute the X-ray temperature from the spectra 


o 

< 



-2 


Arc Minutes 


Fig. 9. — A 4’x4’ temperature map of RCS2327, derived from 
the deep Chandra Cycle 13 data. The image is centered on 
RA=23:27:27.53, Dec=—02:04:35.6. The temperature is indicated 
by the color scale, in keV. The contours are X-ray flux isophotes. 
Point sources have been masked. The temperature and flux maps 
clearly indicate that RCS2327 has a cool core, and a regular mor¬ 
phology, with no significant substructure. 

extracted in the spherical shell with maximum and mini¬ 
mum radii equal to the specific radius and 15% its value. 
The obtained Yx is compared with the Yx — M relation 
calibrated from hydrostatic mass estimates in a nearb y 
sample o f clusters observed with Chandra (Vikhlinin et| 
ah|2009 ). This returns an estimate of M 500 and, thus, a 
new value for i? 5 oo- The process is repeated until con¬ 
vergence in the radius estimate is reached. The result¬ 
ing radius is R 500 = 1.27+ 0 '5g/i7o Mpc, corresponding to 
M 50 o = 1.27t°;® 4 x 10 ls /ifg L M o . 

While the two extrapolation methods agree within er¬ 
rors, we note that the NFW-mass extrapolation results in 
much larger uncertainty. This is due to the fact that the 
HSE mass profile is constrained at a small radius (just 
above R 2500 ) and thus the external slope of the cluster 
is poorly constrained from X-ray observations. 

3.4. Sunyaev Zel’dovich Array Observations and SZ 
Mass Estimates 

The Sunyaev Zel’dovich Array (SZA) observed 
RCS2327 for a total of 48 hours between 2007 Septem¬ 
ber and November. The SZA is an eight-element inter¬ 
ferometer with 30 and 90 GHz receivers. The SZA was 
configured in a standard configuration with 6 telescopes 
arranged in a compact, short-baseline configuration with 
two outlying telescopes ^30 nr from the central group. 
The correlated bandwidth was 8 GHz, centered on 31 
GHz, resulting in projected lengths of 350—1300A on the 
short (SZ-sensitive) baselines and 2000—8000A on the 
longer baselines. The data were calibrated an d flag ged 
using the MATLAB pipeline described in Muchovej et 
ah] (|2007); 41% of the data were removed, largely due 
to shadowing in the compact array, which is increased in 
equatorial and lower declination objects. The rms noise 
level in the short-baseline data is 0.21 nrJy beam -1 , cor¬ 
responding to a 15 /iK rms brightness temperature in 
the 1.9'x 2.8' synthesized beam. A bright radio source 
is detected nearly 9' to the west of the cluster, but it 
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Fig. 10. — Radial profiles derived from the analysis of Chandra Cycle 13 150 ks X-ray observation, extended out to 800 kpc. (a) Top-left: 
Radial profile of the X-ray surface brightness in the soft band [0.5-2.5] keV (data points). The dotted line represents the best fit model, 
(b) Top-right: Projected temperature profile, clearly showing a temperature decrease towards the center of the cluster, indicating a cool 
core. The dotted line represents the best fit model, (c) Bottom-left: Metallicity profile (d) Bottom-right: Hydrostatic mass profile (solid 
line) and its 1 -a uncertainty (shaded area). The vertical line indicates i?2500- 


does not affect the SZ detection. There is one faint 
radio source coincident with the cluster which we jointly 
model when presenting data. This source is present in 
the NRAO VLA Sky Survey (NVSS, Condon 1998). The 
deconvolved image of the cluster, after subtraction of the 
radio sources, is shown in Figure [IT] The peak signifi¬ 
cance in this image is 22 <r. 

The SZA interferometer acts as a spatial filter sensi¬ 
tive to the Fourier transform of sky emission on angu¬ 
lar scales determined by the baseline lengths. To re¬ 
cover the integrated Compton-j/ parameter, Y, we fit our 
Fourier plane data to the transform of t h e gen e ralize d 
NFW pressure profile presented in Nagai et al. (2007), 
which was motivated by si mulations and X-ray observa- 
tions (Mroczkowski et al. 20091. In this five-parameter 
model we fix the three shape parameters ( a , /3, 7 ) to the 
best fit values derived from X-ray o b servations of clus¬ 
ters (1.0620, 5.4807, 1.156; |Arnaud et al.||2010|). We fit 
the profile normalization and scale radius. The cluster 
centroid and the flux of one emissive source are allowed 


to vary as well. The models are fit to the data directly 
in the uu-plane, which correctly accounts for the noise in 
the data. 

To determine the significance of the cluster detection, 
we compare the \ 2 of the best fit model including the 
cluster and emissive source with the \ 2 °f the best fit 
model including only the emissive source. Expressed in 
terms of Gaussian standard deviations, the significance 
of the SZ detection is 30.2o\ 


3.5. Estimates of the Y parameter 

We compute two estimates of the Y parameter, Y sp h 
and Y cy i. Y sp h is a spherical integral of the pressure pro¬ 
file. It is relatively insensitive to unconstrained modes 
in the interferometer-filtered data and proportional to 
the total integrated pressure of the cluster, making it a 
robust observable (see Marrone et al 2011). 

To compute Y sp h , we volume-integrate the radial pro- 
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file to an overdensity radius, va, 

rr a 


I A ,sph — 


<7 T 

m e c 2 


P(r/r s )dV, 


( 8 ) 


/o 


as in Marrone et al. (2012). We determine the overden¬ 
sity radius ot integration by enforcing consis tency with 
the Ysnn.snfo — M^no scaling relation derived by Andersson 
et al. (2011|). To enforce consistency, we iteratively chose 
the integration radius (and, by extension, the mass) un¬ 
til the mass and Ya lie on the mean relation. This 
analysis yields Y 500 ^ sph = 13.4 ± 1.0 x 10 -5 Mpc 2 with 
I? 50 o = 1-13 ± 0.02 Mpc (this radius is 3-15% smaller 
than the R 500 that we derive from extrapolating the X- 
ray data in Section 3.31. 

Y cy i is a cylindrical integral of the pressure profile along 
the line of sight, 


I A ,cyl — 


l7 t 
m e c 2 


dfl 


p e dl 


(9) 


(see also Mroczkowski et al. [2009). Y cy i corresponds to 
the aperture integrated SZ flux, and is sensitive to the 


line of sight contribution of pressure beyond the radius 
of interest. For our gNFW fits, the ratio of Y cy i/Y sp h 
in the 211"8 aperture is 1.0564. We compute Y 500tCy i to 
directly compare our r esults with the SZ obse rvations of 
RCS2327 reported by |Hasselfield et al.|(|2013|) with data 
from the Atacama Cosmology Telescope (ACT). From 
their “Universal Pressure Profile” (UPP) analysis, which 
implicitly imposes the Y — M scaling relation of Arnaud 
(2010), they obtain Y 50 o, cy i = 19.1 ± 0.2 x 10 -5 Mpc 2 
within an aperture of R 500 = 2.8±0.1 arcmin (1.22±0.04 
Mpc at the cluster redshift). Using the same aperture, 
we me asure Y 500 ^ cy i = 1 6.x 10 -5 Mpc 2 , within 1.5er 
of the Hasselfield et al. (2013) measurement. 


3.6. SZ Mass Estimates 

We estimate the cluster mass from the value of Y sp h 
quoted above, which co rrespo n ds to M Rn p = 8. 5 ±1.1 x 
10 14 /iy O 1 Mq using the 


Andersson et al. 


(2011) scaling 


relation. The u ncertainty assumes 21% ~scatter in Y at 


fixed mass in the A nders son et al. (20111 scaling relation 
( Buddendiek et al.||2014 ~ 


We also estimate the cluster mass by applying the 
method outlined in Mroczkowski (2011, 2012) to the SZA 
data. This method assumes the gas is virialized and in 
thermal hydrostatic equilibrium within the cluster grav¬ 
itational potential. Further, we assume the total matter 
den sity p tn t follows an NFW dark matter profile (Navarro 
et al.||1995|), with the gas density p gas is a constant trac- 
tion ot the total density (p ga s = / gas ptot), and the pres¬ 
sure and density profiles are spherically symmetric. This 
method has been applied successfully and co mpare d with 


2012 

Umetsu et al. 

2012 

Medezinski et a 

1 

2013). A 

nt to a gNFW prbfi 

e described above yieJ 

d 

s A 2500 = 


0.54T0.01 Mpc, M 2500 = 4.6±0.5 x lO 14 /if o L M 0 , assum¬ 
ing an average gas fraction wi thin -R 2500 of f gas = 0.137 
from the x-ray analysis in § |3.3| This method yields 
R 500 = 1.15± 0.04 Mpc, M 500 = 8.9 ±0.9 x lO^/i^M©, 
assuming an average gas fraction in this radius of f gas = 
0.12 from Menanteau et al. (2012). We estimate a ~ 10% 
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Fig. 11.— SZ observation of RCS2327. The CLEANed image is 
8 ; x 8 ' in size. The image is in units of SNR and made from the 
SZA uu-data with a Gaussian uu-taper of 0.1 at 4k A. The FWHM 
of the synthesized beam for the image is shown in the bottom left 
corner. In this image, SNR of 1 corresponds to 0.21 mJy/bm. 

scatter due to the assumption on average gas fraction 

value and other model assumptions. _ 

Our mass estimates are consistent with Hasselfield 
et al. (2013), who measure M 50 o 
TUU-i 


= 9.4 ± 1.5 x 


ft 7o M 0 from the UPP Y pa rameter qu oted above. 


10 1 

In addition to the UPP mass, Hasselfield et al. (120131 


report a range of higher M 500 estimates based on differ¬ 
ent scaling relations, M 50 o = 12.5 — 14.3 x lO 14 /iy O 1 M 0 , 
somewhat higher than our measurement. However, the 


incons istency between the higher-mass Hasselfield et al. 
(20131 values and our measurement is not significantly 


worse than the inconsistency with their own UPP mass. 

3.7. Wide Field Imaging and Weak-Lensing Mass 
Estimates 

We obtained deep wide-field imaging data for RCS2327 
in the i! filter using Megacam on CFHT with the aim of 
determining the mass using weak gravitational lensing. 
The observing strategy and weak lensing analysis fol¬ 
lows tha t of the Canadian Cl uster Comparison Project 
(CCCP; |Hoekstra et al. 2012), with the only difference 
that we use the i' tor the weak lensing analysis. The i' 
data consist of 8 exposures of 650 s each, which are com¬ 
bined into two sets (each with a total integration time of 
2600 s). The pointings in each set are taken with small 
offsets, such that we can analyse the data on a chip-by¬ 
chip basis. 

The various steps in the analysis, from object detec¬ 
tion to unbiased shape m easu r eme nts and cluster mass, 
are described i n detail in Hoeks tra (2007), with updated 
procedures in Hoekstra et al. 1 2015|) and we refer the 
reader to those papers for mo re details. We measure 


galaxy shapes as described in Hoekstra et al. (2015), 
which includes a correction for multiplicative bias based 
on simulated images. The resulting shapes are estimated 
to be accurate to 1 — 2%, much smaller than our sta¬ 
tistical uncertainties. The shape measurements for each 
set of exposures are then combined into a master cata- 































































































12 


Sharon et al. 



0 500 1000 


r [”] 


Fig. 12.— Mean tangential shear as a function of radius from the 
BCG is shown in the top panel, along with the best fitting isother¬ 
mal sphere model for reference. The model was fit only to the 
data at r > 75” (solid line). The dashed line is the extrapolation 
of the model to smaller projected radii. The bottom panels shows 
the ’B’-mode lensing signal, which should be consistent with zero 
if systematic distortions have been correctly removed. 


log which is used to derive the weak lensing mass. To 
reduce contamination by cluster members, we also ob¬ 
tained four 720 s exposures in r', which are combined 
into a single image. Galaxies that are located on the 
cluster red-sequence are removed from the object cata¬ 
log, which reduces the level of contamination by a factor 
of two. However, many faint cluster members are blue, 
and we correct the lensing signal for this r esidual con¬ 
tamination, as described in Hoekstra (2007). 

To quantify the lensing signal, we compute the mean 
tangential shear as a function of distance from the clus¬ 


ter center using galaxies with 22 < i' < 24.5. Figure 12 


shows the resulting signal, which indicates that the clus¬ 
ter is clearly detected. The bottom panel shows a mea¬ 
sure of the lensing ‘B’-mode, which is consistent with 
zero, indicating that the various corrections for system¬ 
atic distortions h ave bee n properly applied. 

As discussed in Hoekstra (2007), the weak lensing mass 
can be derived in a number of ways. However, to re¬ 
late the lensing signal to a physical mass requires knowl¬ 
edge of the redshift distribution of the galaxies used i n 
the lensing analy sis. We use the results from [Hoek¬ 


stra et al. (2015) and find that the mean ratio of an¬ 
gular distances between lens-source and observer-source 
is Di s /D s = 0.194. 

For reference, we show t he b est fit singular isother¬ 
mal sphere model in Figure 12 for which we obtain an 
Einstein radius Re = 10 , /1 ± lv 9, which yields a velocity 
dispersion of a = 1345(^34 km s -1 for the adopted source 
redshift distribution. This value is in excellent agreement 
with the dynamics inferred from the galaxy redshifts. We 
also fit an NFW model to the data, adopting the mass- 


concentration relation suggested by Duffy et al. (2008), 
which yields a mass M 2 oo = 2.0) l Io)8 x IO^/i^Mq. We 
compare the weak len sing mass to other mass estimates 
and in other radii in § |4.2[ 

3.8. Strong Lensing Mass Estimates 

RCS2327 was observed by HST- l-ACS (Cycle 15 pro¬ 
gram GO-10846; PI Gladders) as part of a larger effort 
using both ACS and NICMOS to acquire deep multi¬ 
band imaging of this cluster. Unfortunately, the fail¬ 
ure of ACS in early 2007 truncated this program, and 
the only complete image which was acquired is a 3- 
orbit F435W image of the cluster core taken using the 
ACS Wide Field Channel^ Additional available obser¬ 
vations of the cluster core relevant for the strong lens¬ 
ing analysis include a deep (~2 hours) K-band image of 
RCS2327 acquired using the PANIC instrument on the 
Baade Magellan I telescope in 2006, as well as an incom¬ 
plete 4-pointing mosaic of RCS2327 in the F160W filter 
taken with HS T+ NICMOS; we have reconstructed this 
last image from the useable portions of a nominally failed 
HST observation which nevertheless yielded some useful 
frames in a single orbit before guiding issues truncated 
the remainder of the observations. A color composite 
image of the cluster core, made fro m the F435W image, 
the deep LDSS-3 j-band image (see § 3T] above), and the 
PANIC K s -band image, is shown in Figure [13] 

Using these various imaging data, we identify two sets 
of multiply-imaged galaxies that are lensed by RCS2327 
for which we have acquired spectroscopic redshifts as 
p art of the overall spectroscopic program des cribed in 
§ |3.2[ Both sources are indicated in Figure |7TT| A merg¬ 
ing pair of images of source A is located at 23:27:29.43, 
—02:03:47.8, to the NE of the cluster center. Its redshift, 
z = 2.9834 ±0.0010, is determined from a strong Lya 
emission line present in the early LDSS-3 spectroscopy 
described above. This lensed source was apparent in the 
RCS2 discovery imaging data, with a remarkably large 
separation from the cluster center, R=56"8, as measured 
from the BCG. The arc does not appear to be caused by 
local substructure in the cluster, as there are no nearby 
significant cluster galaxies. 

Source B was observed spectroscopically in queue mode 
in semester 2007B using the Gemini South telescope 
with the GMOS instrument. We observed RCS2327 for 
8 x 1800s in multi-object spectroscopy mode. The obser¬ 
vations were taken with the B600_G2353 grating, no fil¬ 
ter, and the detector binned 1 x 2 (spatial x spectral axes), 
resulting in wavelength coverage of ~ 2700A per slit, 
and a spectral resolution of ~ 240 km sec -1 . The grat¬ 
ing tilt was optimized to record a wavelength range of 
~ 3800 — 6500A for images of source B. 

The redshift of source B is z = 1.4155 ± 0.0008, based 
on [O IIJAA3727 emission present in a GISMO observa¬ 
tion (see § 3.2) confirmed by several Fell lines in absorp¬ 
tion in the Gemini spectra. Source B is lensed into three 
images, and is not morphologically obvious in the discov¬ 
ery data from RCS2327, as it is not lensed into a classic 
tangential arc. It is apparent in the combined HST and 


22 The field of RCS2327 was recently imaged by HST in Cycle 
20 program GO 13177 (PI Bradac). These data are not used in 
this paper. A forthcoming lensing analysis of the Cycle 20 data 
will be presented in|Hoag et al.|(2015|. 
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Fig. 13.— A 2’x2’ color composite image of the core of RCS2327 composed of images from HST in the F435W filter (blue channel), 
and from Magellan in the i —band (green channel) and ff s -band (red channel). The point spread functions have been matched to the worst 
image; the effective resolution is ~ 0'.'6. The merging pair image of source A at z= 2.9834 is indicated. The greyscale cutouts are 2" in 
diameter and show the full-resolution images of source B, at 2=1.4155, in the F435W filter. The critical curves from the best fit lens model 
are overplotted, in red for a source at 2=1.415, and in orange for a source at 2 =2.9834. 


IR imaging since it has a unique color and internal mor¬ 
phology. These properties also allow us to robustly elim¬ 
inate the presence of a fourth counter image; source B 


is lensed as a naked cusp configuration (e.g., Oguri & 
Keeton 2004). A close inspection of the F435W image 


reveals that two of the images (B1 and B2) have two 
emission knots at their center; overall source B appears 
to be a compact galaxy with a primarily redder stellar 
population, but with two well confined regions of active 
star formation in the galaxy’s core. The detailed position 
of these bright knots indicates a larger magnification in 
the tangential direction than in the radial direction for 
this source. The two knots in the third image are not 


resolved, but the image is elongated in the tangential di¬ 
rection. Source B also has a significant Einstein radius, 
with separations for the three images from the cluster 
BCG of 3678, 3676, and 3578. Further lensed features 
are also apparent, but we do not yet have redshift in¬ 
formation for them and they are not used in the initial 
lensing model discussed below. 

A strong lensing model for RCS2327 was constr ucted 
usin g the_publicly available software LENSTOOL (Julio 


The mass model is composed of muiti- 
umps. The cluster h alo is represented by 


et al. 2007). 
pie mass c 

a generalized NFW distribution (Navarro et al. 1997) 
parametrized with position, a;, y\ ellipticity e; position 
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angle 6 ; central slope a; and concentration c. The 50 
brightest red-sequence-selected cluster-member galaxies 
are represented by Pseudo- Isot hermal Ellip soidal Mass 
Distributions (PIEMD; see Julio et al.|2007 for details) 
parametrized with positional parameters (x, y , e, 9) that 
follow their observed measurements, r core fixed at 0.15 
pc, and r cu t and cto scaled with their luminosity (see 
Limousin et al. 2005 for a description of the scaling re¬ 
lations). The parameters of an L* galaxies were fixed at 
Tout = 40 kpc and <7o = 160 km s -1 . The model consists 
of 13 free parameters. All the parameters of the cluster 
halo are allowed to vary (R.A., Deck of the mass clump, 
ellipticity, position angle, scale radius, concentration and 
central radial mass profile). 

The constraints are the positions of the lensed features 
and their redshifts. Each component of Arc A was repre¬ 
sented by three positions, and the two cores of source B 
were used in each of its images. The best fit model is de¬ 
termined through Monte Carlo Markov Chain (MCMC) 
analysis through minimization in the source plane, with 
a resulting image-plane RMS of O'/17. The best fit pa¬ 
rameters and their 68% percentile uncertainties are pre¬ 
sented in Table IT] Some of the model parameters are not 
well-constraincaby the lensing evidence. In particular, 
a large range of values is allowed for r s and a, and the 
model can converge on any value of the concentration 
parameter c. The latter is not surprising, since in order 
to determine the concentration parameter one needs to 
constrain the slope of the mass profile on small and large 
radii, beyond the range of the strong lensing constraints. 
Thus the concentration uncertainty given in Table [l] rep¬ 
resents the range of priors assumed in the lens modeling 
process. We find strong correlations between, a, r si and 
c, which we fit to find r s = 575.4 —906.8a+ 583.3a 2 and 
c = 16.1 — 8.8a. 

The Einstein radius of a lens is often used as a mea¬ 
sure of its lensing cross section, or strength. We mea¬ 
sure the effective Einstein radius as Re = \JA/'k 1 where 
A is the area enclosed by the tangential critical curve, 
Re{z = 1.4155) = 25'.'9 for source B, and Re{z = 
2.9834) = 40"2 for the giant arc A. These radii are 
smaller than the separations between the arcs and the 
BCG, due to the ellipticity of the lensing potential. The 
mass that results from the lensing model can be quoted 
at a range of radii, though it is clear that the mass 
is most robustly measured at the critical radii probed 
by the lensed images us ed to construct the model (e.g. 
Meneghetti et al. 2010). We integrate the strong lens 
model within circular apertures at radii corresponding 
to the mean positions of sources A and B with respect 
to the position of the main cluster NFW halo, and find 
enclosed projected mass M cy i(< Ra) =5.8t(j;2 x 10 14 M o , 
M cy i{< Rb) =3.4+qj x 1O 14 M 0 . Statistical uncertain¬ 
ties are computed by sampling models described by the 
MCMC outputs, considering only models with values of 
X 2 within two of the best fit, representing 1-cr uncer¬ 
tainty in the parameter space. The resulting masses are a 
measure of the projected (i.e., cylindrical) masses within 
the quoted radii. These statistical uncertainties may fail 
to reflect some systematics due to the small number of 
lensing constraints in this system. In particular, since 
the lens is only constrained by arcs on one side of the 
cluster, we see correlations in the parameter space be¬ 


tween the mass, ellipticity, and position of the lens. The 
superior data expected from HST Cycle 20 program GO- 


13177, will enable a better constrained lens model (Hoag 
et al. 2015). We adopt a 15% systematic uncertainty 
from Zitrin et al. (2015) for clusters with similar strong 
lensing signal. 

A notable further result from the strong lens model 
is that the cluster halo is offset from the BCG by 
1-17-CL24 an d 7.42tcK63 arcseconds in right ascension 
and declination, respectively. This corresponds to an 
offset of 54 kpc at the cluster redshift. Figure [5] shows 
the positional relationship between the cluster galaxies 
- as demarcated by red-sequence members - and both 
the X-ray data and the strong lensing model. The peak 
of the X-ray emission is coincident with the position of 
the BC G as is typically seen in lower redshift re laxed 
clusters (Sanderson et al. 12009 Bildfell et al.|2008|. The 
center of the overall distribution of the red sequence light 
is coincident with the strong lensing mass peak, both of 
which are hence offset from the BCG and the X-ray cen¬ 
troid by ^60 kpc. Disagreements between the mass peak 
as traced by lensing and th e X-ray centroid are seen in 
major clusters mergers (e.g. Bradac et al.|2008 Mahdavi| 
et al|2007||Clowe et al.|2004p although the magnitude of 
the disagreement in RCS2327 is not nearly as large and 
is similar to that observed in inter mediat e cooli ng flow 
clusters in the sample of clusters in Allen (1998). How¬ 
ever, the differing positions indicated by various mass 
tracers is arguably the strongest evidence that RCS2327 
is anything but a single relaxed halo; we explore the im¬ 
plications of this further in § [4j 

4. DISCUSSION 

The various mass proxies detailed in § [3] all indicate 
that RCS2327 is an exceptionally massive cluster given 
its redshift. Each of these mass proxies is naturally sen¬ 
sitive to the mass of RCS2327 at a particular radius, 
and involves in all cases one or more simplifying assump¬ 
tions that allow the conversion of the observable signal 
into a mass estimate. For example, the X-ray data most 
directly constrain the mass at an overdensity radius of 
an< J conversion of the X-ray spectrum and radial 
luminosity profile to a spherical mass estimate requires 
the assumption of hydrostatic equilibrium. The strong 
lensing data are sensitive to mass at similar or smaller 
radii than the X-ray data, but fundamentally measures a 
cylindrical mass in projection. The galaxy dynamics are 
sensitive to mass at the virial scale and rely on external 
scaling relati ons to provide a mass estimate, which, as 
detailed in § 3.2 are sensitive to not well known issues 
of velocity bias and orbital anisotropies. 

4.1. Comparison to Other Clusters 

The left and middle panels of Figure [14] compare the 
velocity dispersion, X-ray temperature and richness of 
RCS2327 to the global correlations of these properties in 


an intermediate X- ray selected cluster sample from Yee 
& Ellingson (20031. The right panel of Figure 14 plots 


the measured velocity dispersion and X-ray temperature 
of R CS232 7 against th e cluster data and fitted relation 
from Xue & Wu (2000). We plot both the total and red- 
sequence richnesses, and the velocity dispersion from all 
galaxies, and only early-type galaxies. Which of each of 
these properties is best compared to the correlations in 
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TABLE 1 

Best fit Strong Lensing Model Parameters 


Halo 

Model 

RA 

Dec. 

e 

0 

r s 

a 

C 



(") 

(") 


(deg) 

(kpc) 



Halo 1 

gNFW 

1 17+0.47 

-0.24 

7 42~*” 1,42 

0 26 - *” 0 ’ 04 
u - zd -0.06 

1.02. 

486 ±“ 6 

1 44+0.09 
■ L *”— 0.62 

Q 4Q+5.09 
o.^U—o 40 


Note. — Coordinates are measured in arcseconds East and North of the center of the BCG, at [RA, Dec]=[351.865026, —2.076924], 
The ellipticity is expressed as e = ( a 2 — ft 2 )/(a 2 + b 2 ). 6 is measured North of West. Error bars correspond to 1 -a confidence level as 
inferred from the MCMC optimization. 



Fig. 14.— The comparison of richness, velocity dispersion, and X-ray temperature for RCS2327 (green s quares') against re l ations for 
these quantities from the literature. The left and middle panels compare to the data and fitted relations in|Yee & Ellingson| (2003) for 
ve locity d isper sion - richness (left) and X-ray temperature - richness (middle) and the right panel compares to tne data ancTfitted relation 
in Xue Wu (2000) for velocity dispersion - X-ray temperature. Both the total and early-type-only richnesses are shown, as well as 
velocity dispersions from both the early-type-only and all galaxies. The X-ray temperature and velocity dispersions (for all galaxies and 
early-type-only) of ACT-CL J0102-4915 (z = 0.87) from[Menanteau et al. i|2013[) are plotted in orange circles. 


Yee & Ellingson (1200 3b or Xue & Wu (20001 is not ob- 
vious (e.g., see |3.1| ). Regardless, to within both the 
measurement uncertainties and these systematic uncer¬ 
tainties, these three measures (which probe large scale 
dynamics, the gas properties of the cluster core - and 
hence small scale dynamics, the gas fraction, and the 
like - and the stellar mass-to-light ratio) are all consis¬ 
tent with a massive cluster with properties drawn from 
the global correlations seen in lar ge c l uster sampl es. 

In 


Figure [l5| reproduced from 


plot the estimated M 50 o versus the redshitt o: 


Bleem et al. 


(2015), we 
RCS2327, 


compared to clusters from large X-ray and SZ cluster sur¬ 
veys. The figure illustrates that RCS2327 is among the 
most massive clusters at all redshifts, and in particular 
at z > 0.7. 


4.2. Comparison of Mass Proxies 

Though each of the mass proxies discussed above mea¬ 
sure the mass most naturally at differing radii, it is still 
instructive to compare the results directly. To do so we 
consider several additions to the main analyses in § [3j 
Table ED summarizes the mass estimates from different 
mass proxies at different radii, and they are plotted in 
Figure [TO) 


SPT-SZ 2500 deg 2 X 
ROSAT-AII sky ■ 
Planck-DRI ♦ 
ACT 

RCS2327A 



Fig. 15.— Reproduced from Bleem et al. |2015|). Estimated 
mass versus redshift for cl usters fro m lour large X-ray and SZ sur¬ 
veys :^_SPT2SZ_2500_deg2_J Bleem et al.bo iHt, ROSAT all sky sur¬ 
vey dPiff aretti et al.|| 20lT , Flanck-UKl ([Pl anck Collaboration et| 
|al.||2UT3| ), and AUT (Marriage et al. 2 011} . The X-ray mass ol 
RCS2327 (M 500 — 1-3^0 2 x 10 15 M©) is overplotted as a ma¬ 

genta triangle, placing it among the most massive clusters across 
all redshifts, and comparable to only few other clusters at z > 0.7. 
We note that clusters that appear in different catalog may show 
multiple tim e s on this plot, e.g., ACT- CL J0102—4915 (|Marriage| 
|et al.| (|2011|) , |Menanteau et al.| (|2013|) ) appears with three mass 
estimates at z = U.87. Utiier notable nigh-mass clusters are SPT- 
CLJ2337-5942 at 2 = 0.77, SPT- CLJ0615-5746 at z = 0.972 and 
SPT-CLJ2106-5844 at z = 1.132 iBlee m et al.ll2015l |Foley et al.| 
|2011| |Planck Collaboration et al.|20fT| |. 


4.2.1. Cylindrical Masses from X-ray, Strong and Weak 

Lensing 

As noted above, weak and strong lensing are both sen¬ 
sitive to projected mass density. However, they probe 
different regimes of the mass distribution: strong lensing 


is insensitive to the mass at the outskirts of the clus¬ 
ter, where no strong-lensing evidence exists. Weak lens¬ 
ing lacks the resolution at the cluster core. To com¬ 
pare the weak and strong lensing mass estimates, we 
first compute the projected enclosed mass (also known 
as the aperture mass) as a function of radius directly 
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fro m the weak lensing data. We u se the Cc statis¬ 


tic (Clowe et al. 


1998 Hoekstra 2007) and convert the 


measurements into projected masses, using the best fit 
NF W to estimate t he large scale mean surface density 
(see|Hoekstra|[20071 for details). The dependence of the 
resulting projected m ass estimate on the a ssumed den¬ 
sity profile is minimal (Hoekstra et al.|2015). At a radius 
of 500 kpc this yields a projected enclosed mass of 
M W L,cyi(< 500/ifokpc) = 5.7 ± 1.1 x lO^/i^Mg. 

We can similarly extend the mass estimate from the 
strong lensing model to larger radii. However, since the 
lens model is only constrained by lensing evidence in the 
innermost 400 kpc (measured from the BCG) we in¬ 
crease the systematic uncertainty of the strong lensing 
mass esti mate by ~10%. Following the analysis out¬ 
lined in § 3.1 we find a mass at a SOOfi^/kpc radius of 
A/sL, C yi(< SOO/i^Q 1 kpc) = 8.0 ± 1.5 x 10 14 ft- 7 O 1 M Q . These 
two values are in fair agreement. We refrain from extrap¬ 
olating the strong lensing mass to larger radii, where the 
strong lensing model is not constrained. 

The X-ray masses can be converted to cylindrical mass, 
by integrating along the line of sight out to 10 Mpc on 
both sides of the cluster center. We note that this may 
introduce some uncertainty as this is model-dependent. 
The projected enclosed X-ray masses at the radii of 
the lensed galaxies (see Table [ 2 J) are 2.01 q;|, 4.0lJ;g, 
6 . 2 ti; 4 xl 0 14 /i 7 O 1 M Q for 271, 352, 500 kpc, respectively. 
These values are in fair agreement with the projected 
enclosed masses derived from strong lensing, 3.4j(Q'' 4 , 

5. 8 ^ 0 ( 2 ; 8 . 0 (tg; 4 xl 0 14 /i 7 O 1 MQ, respectively. The differ¬ 
ences are in line with expected uncertainties and biases 
(see, e.g., Mahdavi et al. 2013) for hydrostatic masses, 
as overall we find that the lensing masses are somewhat 
higher than the X-ray and SZ masses. Nevertheless, it 
may also indicate that structure along the line of sight or 
elongation of the cluster halo may be significant. For ex¬ 
ample, the structure that is indicated by a concentration 
of galaxies at z ~ 0.73 (Figure 7] may be contributing 


ture lensing analysis ( 

D’Aloisio et al. 

2014 

McCully et 

al.|2014 

Bayliss et al. 

20 MaF 


4.2.2. Spherical Masses 

To compare the weak lensing, X-ray, and S Z masses 
we de project the aperture masses following IHoekstra 


d2007), assuming the mass-concentration from Duffy et 


al. (2008). Although the deprojection is somewhat model 


dependent, it is less sensitive to deviations from the NFW 
profile. At the cluster core, we compute the correspond¬ 
ing deprojected weak lensing mass within 500 hfij kpc 
(approximately i? 2 soo)- We obtain a value of l(< 
SOO/i^kpc) = 4.1+j'^ x 10 14 /i 7 0 1 Mq within this ra¬ 
dius, in agreement with the X-ray estimate of Mx ,2500 = 
3. 2 ^ 0)3 x 10 14 /i7 0 1 Mq, and SZ mass of Msz,2500 = 4.6 ± 
0.5 x 10 14 /i to 'Mq. 

At large radii, we u se the extrapolated X-ray mass 
as described in § |3.3| In making this comparison we 
note that the native values of i? 2 oo from each of these 
analyses agree within the uncertainties. The X-ray mass 
is Mx, 2 oo = l- 8 (to .7 x 10 15 /ito m O an d the weak lens¬ 
ing mass from the NFW fit is M\v l, 2 qo(< 2.1Mpc) = 


2.7±0.7 x 10 15 /i7 O 1 Mg. Hence at large radii the extrapo¬ 
lated X-ray mass and weak lensing data also agree within 
the uncertainties. 


4.3. Mass Profile 

Figure [16] presents the enclosed masses measured in 
this paper as a function of cluster-centric radius, as 
well as SZ masses from the literature. As demon¬ 
strated above, these measurements are consistent with 
each other within errors, and trace the mass profile from 
the very core out to i? 200 - 

We fit a set of spherical NFW profiles (Eq. [7]) to the 
spherical masses measured in this paper. To estimate 
the range of fits that are consistent with the measure¬ 
ments, we fit the profile 1000 times, each time to a set 
of measurements that were randomly sampled from their 
l-u uncertainties, and weighted by their uncertainties. 
We did not include in the fit the extrapolated estimates 
and constraints from the literature. These masses are 
shown in Figure [l6| for reference, extrapolated measure¬ 
ments with in dashed error bars, and the Hasselfield et 
al. (2013) mass estimates in thin lines. A large range 
of scale radii is consistent with the measured masses, 
and the resulting range of NFW profiles is shown as the 
solid shaded area in Figure [l6| The striped area in the 
right panel of Figure [16] traces the cylindrical mass from 
the same NFW profiles that were fit to the spherical 
masses. While we could simultaneously fit the profile to 
the cylindrical and spherical masses, we choose not to, 
because the cylindrical strong lensing measurements do 
not assume spherical symmetry and thus should not be 
expected to be described by a spherical NFW profile. We 
find that the strong lensing masses are somewhat higher 
than the predicted cylindrical masses, which could be due 
to the triaxiality that is not taken into account in this 
simplified fit. As expected, the projected X-ray masses 
do agree with the spherical profile, since they were com¬ 
puted by integration of the X-ray best fit spherical profile 
along the line of sight. 

The simplistic NFW fit to all the non-extrapolated 
cylindrical mass measurements yields r s = O^iJ ^Mpc. 
However, while a fit of a spherical NFW profile to the 
mass measurements is possible (though a large range of 
scale radii is consistent with the results), we caution that 
such a fit is not meaningful at this point. The different 
measurements were conducted completely independently 
of each other, and rely on different assumptions as de¬ 
scribed in the previous sections (e.g., mass-concentration 
relations, spherical symmetry, hydrostatic equilibrium, 
various scaling relations). In particular, some of the 
mass proxies already assume a certain mass profile slope. 
A self-consistent combined multi-wavelength analysis is 
called for. Such an analysis would ideally allow triaxial 
symmetry, and fit the mass distribution simultaneously 
to constraints derived directly from all the observables: 
strong lensing constraints, weak lensing shear, galaxy 
velocity distribution, and X-ray and SZ measurements. 
This sort of analysis is left for future work, and is not 
within the scope of this paper. 

5. SUMMARY AND CONCLUSIONS 

We present a multi-wavelength analysis of RCS2327, 
a massive cluster at z = 0.6986. The mass is estimated 
independently at several radii, using seven different mass 
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TABLE 2 
Estimated Masses 


Mass proxy 

Projected Mass [10 14 

hjQ Mq] 



Spherical Mass t 




R = 217kpc 

R = 352kpc 

R = 500kpc 

T2500 

A4250o/10 14 

T500 

M 500 /IO 14 

P200 

M 20 o/ 10 14 





[kpc] 

[h ro Mg] 

[Mpc] 

[h 70 Mq] 

[Mpc] 

[h 7 Q Mq] 

Strong Lensing 

3.4l°f±0.5 

5.8±°' 4 ±0.9 

[8.0l°®±1.4] 







X-ray 

2 q+ 0.5 
z ' u -0.4 

4 0 +1 ' 2 
^- U - 0.8 

u - -1.4 

4711|| 

q Q + 0.6 

°- z -0.3 


[lit?] 

[l-78to.43l 

[18±^ s ] 

Weak Lensing 



5.7 ±1.1 

517 

4.3 ± 1.2 

1.34 

ist?;? 


20 t® 

SZ 






1.13 ±0.02 

8.5 ±0.4 


SZ (Mil*) 

Velocity dispersion 




540 ±6 

4 6” 1 ” 0 ' 1 

1.15 ±0.02 

8.9 ±0.8 

2.1 

29-7 ±& 

Caustics 








2.1 

29±7° 

Richness 








2.1 

321^ 


Note. — Summary of the mass estimates from the different mass proxies considered in this work. Square brackets indicate extrapolated 
values. Projected X-ray mass was computed by integrating the mass model along the line of sight out to 10 Mpc on both sides of the 
cluster. 

T The different mass proxies were estimated within different radii, as indicated. 

^ SZ measurments using the method of Mroczkowski (2011). 



r [kpc] R [kpc] 


Fig. 16.— The mass estimates from the different mass proxies considered in this work (see Table |2|l are plotted as a function of radius, 
color-coded by mass prox y as in di cate d in the legend. Points with dashed error bars are from extrapolated results (see text). The SZ mass 
estimates from Hasselfield et al.|(2013| are plotted at R = 1.2 Mpc in green squares. In the left panel we plot spherical masses within 
radius r, and in tne ngnt panel are cylindrical (projected) masses enclosed within projected radius R. 1 -rr uncertainties are shown; we 
note that when and M& are determined jointly their uncertainties are correlated. The shaded area in the left panel is the 1 -rr range 
of spherical NFW mass profiles that were fit to the spherical masses measured in this work. The measurements that were included in the 
fit are indicated with thick circles and errorbars. The cylindrical masses in the right panel were not include d in the fit, nevertheless, we 
show the projected mass density of the same fits as striped area in the right panel. As discussed in Section |4.3[ this simple fit does not 
represent a true joint analysis of the data, since various assumptions on the slope of the mass profile are already folded into some of these 
measurements. 


proxies. At the core of the cluster, we measure the pro¬ 
jected mass from a strong lensing model; at intermedi¬ 
ate radii, ~ 0.5 Mpc, the mass is estimated from X- 
ray, weak lensing, and the SZ effect. At large radii, 
~ 2 Mpc, we measure the cluster mass from its weak 
lensing signal, the dynamics of galaxies in the cluster, 
and from scaling relations with the richness of the clus¬ 
ter. This analysis provides a unique opportunity of com¬ 
paring methods and testing them against each other at 
a significant redshift. In the previous section we com¬ 
pared mass estimates at overlapping radii. Each of the 
mass proxies is prone to statistical and systematical un¬ 
certainties. Moreover, since all the measurements were 
conducted independently from each other, some of the 
mass proxies rely on assumptions (e.g., assumed mass- 


concentration relation or the derived value of r a) that 
are not necessarily uniform among these proxies. This 
unavoidably contributes to the scatter among the derived 
masse s. Nevertheless, the simple internal comparisons in 
§ |4.2[ an d the comparison to global cluster correlations 
in § |4.1| suggests that RCS2327 is not a peculiar object 
(apart from its overall mass) and we thus expect that a 
self consistent analysis would yield results comparable to 
those presented here. 

In summary, all the evidence point to the conclusion 
that RCS2327 is one of the most massive high redshift 
clusters known to date at z > 0.7. 

The set of measurements presented in this paper is 
expected to be improved upon in the near future, with 
deep HST observations that have already been exe- 
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cuted. Further observations will provide constraints for a 
self-consistent modelin g of the three-dimensional cluster 


mass distribution (e.g. 

Umetsu et al.12015 

Limousin et 

al. 2013 Sereno et al. 

1013j), that takes into account the 


effects of triaxiality and orientation on the mass observ¬ 
ables. 
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